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ABSTRACT 

By performing magnetohydrodynamical (MHD) simulations, we investigate the mass loss of 
intermediate- and low-mass stars from main sequence (MS) to red giant branch (RGB) phases. 
Alfven waves, which are excited by the surface convections travel outwardly and dissipate by nonlin- 
ear processes to accelerate and heat the stellar winds. We dynamically treat these processes in open 
magnetic field regions from the photospheres to ~25 stellar radii. When the stars evolve to slightly 
blueward positions of the dividing line (Linsky & Haisch), the steady hot corona with temperature, 
T Ri 10^ K, suddenly disappears. Instead, many hot (~ 10^ K) and warm (> 10^ K) bubbles are 
formed in cool (T < 2 x 10'' K) chromospheric winds because of thermal instability; the red giant 
wind is not a steady stream but structured outflow. As a result, the mass loss rates, M, largely vary 
in time by 3-4 orders or magnitude in the RGB stars. Supported by magnetic pressure, the density 
of hot bubbles can be kept low to reduce the radiative cooling and to maintain the high temperature 
long time. Even in the stars redward of the dividing line, hot bubbles intermittently exist, and they 
can be sources of ultraviolet /soft X-ray emissions from hybrid stars. Nearly static regions are formed 
above the photospheres of the RGB stars, and the stellar winds are effectively accelerated from several 
stellar radii. Then, the wind velocity is much smaller than the surface escape speed, because it is 
regulated by the slower escape speed at that location. We finally derive an equation that determines 
M from the energetics of the simulated wave-driven winds in a forward manner. The relation explains 
M from MS to RGB, and it can play a complementary role to the Reimers' formula, which is mainly 
for more luminous gaints. 

Subject headings: magnetic fields - stars: coronae - stars: late type - stars: mass loss ~ stars: winds, 
outflows - waves 



1. INTRODUCTION 

Solar-type stars have X-ray emitting coronae with the 
temperatures, T > 10^ K, and hot stellar winds stream 
out with the asymptotic terminal velocities, v^o = 100 — 
1000 km s""', comparable to the escape velocities at the 
surfaces. The total mass lost by the wind is negligi- 
ble during the main sequence (MS) phase in compari- 
son with the stellar mass, whereas the wind may play 
a role in the loss of the angular momentunis. Along 
with the stellar evolution to red giant branch (RGB), the 
mass loss rate, M's, increas es by 4-10 orders of magni- 
tude pudge fc Stencel Ill99ll ). while v^o drops to 10-100 
km s""', which is smaller than the surface escape veloc- 
ity. In addition, the temperature of the wind also de- 
creases, which is probably drastic across the "dividing 
Hne" (DL hereafter; Linsky & Haisch 1979) discussed 
below. The mass loss during the RGB phase becomes 
more important, and it may control the later evolution 
of a star itself. For instance, some stars evolve to a pe- 
culiar group of blue horizontal branch stars, and it is 
inferred that it is related to the mass lo ss during the 
RGB phase ()Yi. Demarque. fc Kim1ll997D . In spite of 
its importance, however, the acceleration mechanism of 
red giant winds is not well understood. 

It was reported that cool single RGB stars redward 
of the DL (near spectral type Kl HI in Hertzsprung- 
Russel (HR)-diagram; Figure [1]) showed an evidence 
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Linskv fc Haisch I Il979f) nor soft X-ray (T ^ 
Avres et all 119811 ). IStmcel fc Mullan I (Il980h 
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further pointed out that the stellar winds appear 
cool and massive after stars cross the DL. However, 
iHartmann. Dupree. fc Ravmond I (|1980[ ) observed late- 
type super-giants in the noncoronal zone that exhibit 
both UV emission and blue-shifted chromospheric ab- 
sorption features, indicating that hot gas and cool out- 
flow coexist. Such hybrid stars were also detected by 
later observations (e.g. Ayres et al. 1998), whereas the 
activity levels of UV/soft X-ray are much lower than in 
the coronal zone. These observations show that the sim- 
ple DL picture roughly holds while the reality seems a 
bit more complicated. 

RGB stars are generally slow rotators. The radiation 
luminosities of intermediate- and low-mass stars near the 
DL are not still large, either. Therefore, the centrifugal 
force and radiation pressure are not sufficient in accel- 
erating the stellar winds (e.g. Judge & Stencel 1991). 
Instead, the surface convective layer is expected to be 
the main origin of driving the winds and UV/X-ray ac- 
tivities. Magnetic fields are probably generated by tur- 
bulent dynamo in the convection zone. The energy of 
the surface turbulences is lifted up through the magnetic 
fields, and the dissipation of the energy in upper regions 
leads to the heating of the atmosphere and the accelera- 
tion of the wind. This picture is essentially the same as 
what takes place in the solar corona and wind. Difference 
from the sun is that the surface gravity is much smaller 
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in RGB stars. In a nai've way of thinking, the transition 
from the steady coronae to the cool winds (with some 
hybrid activities) is possibly explained by this gravity ef- 
fect even though other conditions are similar; the cool 
atmosphere comes to flow out before heated up to coro- 
nal temperature (T ^ 10^ K), as the gravity confinement 
weakens. 

Many pioneering works have tried to reveal nature of 
red giant winds (e.g. Hartmann & MacGregor 1980; 
Holzer, Fla, & Leer 1983; Charbboneau & MacGregor 
1995). However, the consecutive processes from the sur- 
face to the wind explained above have not been stud- 
ied self-consistently; previous model calculations started 
from locations far above the photosphere in order to 
avoid difhculties arising from the huge density contrast. 
Moreover, they assume mechanical or phenomenological 
"wave" energies with somewhat ad hoc prescription of 
the dissipation. 

In this paper, we perform dynamical simulations of 
stellar winds from the photosphere to ~ 25 i?* in open 
magnetic field regions, where R^, is stellar radius. We 
give perturbations at the photosphere originating from 
the surface convection. These fluctuations excite waves 
that propagate upwardly. Among various modes the 
Alfven wave travels to a sufficiently large distance to 
contribute to the acceleration and heating of the stel- 
lar wind. We simulate the propagation and dissipa- 
tion of the waves by solving nonlinear MHD equations, 
and treat the heating and acceleration of the gas in a 
fully self-consistent manner without conventional heat- 
ing functions (' ij3.4.2[) . Wind properties, such as M and 
Voc, are determined from the surface conditions with- 
out ad hoc parameters. Another advantage of dynamical 
simulations is that we do not have to care about the 
critical (sonic) point (s), which is often problematic when 
constructing a transonic wind solution under the steady- 
state condition. A stable wind structure is automatically 
selected whether it is transonic or subsonic. 

We alread y carried out simula t ions i n the s olar corona 
and wind (Suzuki fc Inutsuka I l2005l . |2006l, hereafter 
SI05; SI06). We showed that the coronal heating and 
solar wind acceleration in the open field regions are nat- 
ural consequences of the photospheric perturbations; the 
Alfven waves generated from the surface effectively dissi- 
pate in the corona by nonlinear mechanism to sufficiently 
heat and accelerate the solar wind in spite of the energy 
loss via radiative cooling, thermal conduction, and wave 
reflection. The aim of this paper is to study the evo- 
lution of the stellar winds from MS to RGB stages by 
extending this work to cool RGB stars. We investigate 
the stellar winds of a IMq star from MS {i.e. present 
sun) to RGB across the DL, where Mq is the solar mass. 
For comparison we also study a more massive RGB star 
with 3Mq both blueward and redward of the DL. 

Sj2] presents our simulation method. In ^we show the 
results of the simulations, focusing on the evolution and 
time-dependent nature of the stellar winds. We discuss 
comparisons with related works as well as the uncertain- 
ties and limitations of our simulations in 21 

2. SIMULATION SET-UP 
2.1. Simulated Stars 




Fig. 1. — Simulated stars in HR-diagram. The six red stars 
correspond to Models tabulated in Table [T] The black solid lines 
are the evolutionary paths of 1 A/0 and 3 Mq stars with the solar 
abundance (Girardi et al. 2000). The red line with the cross hatch 
is the DL (Linsky & Haisch 1979). The green solid lines indicate 
constant stellar radii, ij* = 1 Rq, 3.1 Rq, 10 Rq, and 31 Rq 
from below. The blue lines indicate constant gravity, logg = 4.4, 
3.4, 2.4, 1.4 from below. The black dashed lines indicate constant 
escape speeds from the stellar surfaces, Wesc.o = 617, 347, 195, 110 
km s~^ from below. The dark green shaded region denotes the 
location of main sequence stars. 

We simulate the Alfven wave-driven stellar winds of 
the various stars summarized in Table [1] and in Figure [TJ 
First, we consider the 1 Mq stars with different stellar 
radii, i?^ = 1 Rq [i.e. the present sun), 3.1 Rq, 10 Rq, 
and 31 Rq, where Rq is solar radius, to study the evo- 
lution of the winds; note that the DL is between the 
latter two cases. We also simulate the stars with 3Mq 
near the DL to compare with the 1 Mq c ases. Here we 
adopt the stellar evolutionary tracks from I Girardi et aTl 
(|200Cl( ) to determine the basic stellar parameters, such as 
and effective temperatures, T^g, which are used to 
estimate the densities, po, at the photospheres and the 
photospheric amplitudes, Svq, of the surface convection 

2.2. Simulation Details 

We consider ID flux tubes that are super-radially open. 
As shown in Figure [21 a certain fraction of the stel- 
lar surface is covered by closed loops so that open flux 
tubes diverge more ra pidly than the radial expansion 
(|KoDD fc Holzer I Il976l ). 

We inject velocity and magnetic field fluctuations orig- 
inating from the surface convection at the photospheres, 
which are the inner boundaries of the simulations ( §2.3p . 
These fluctuations excite various waves, a fraction of 
which propagate upwardly. We treat both energy and 
momentum transfers in stellar winds, namely the propa- 
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Model 






Teff{K) 


log 3 


Dosc,o{km s"^) 


Po(g cm"^) 


(5t)o>(kms-i) 


T-min(min) 


rmax (min) 




R2{R.) 


I 


1 


1 


5800 


4.4 


617 


1.0 X lO-'^ 


1.0 


1 


32 


1.01 


1.2 


II 


1 


3.1 


4800 


3.4 


347 


5 X 10~* 


1.5 


9 


290 


1.031 


1.4 


III 


1 


10 


4500 


2.4 


195 


1.5 X 10"** 


2.8 


87 


2800 


1.1 


1.5 


IV 


1 


31 


3900 


1.4 


110 


4.5 X 10-» 


5.3 


800 


25600 


1.31 


1.8 


V 


3 


17.3 


4700 


2.4 


255 


1.5 X 10"* 


3.7 


90 


2900 


1.1 


1.5 


VI 


3 


54.8 


4200 


1.4 


144 


4 X 10-3 


8.0 


930 


29900 


1.31 


1.8 



TABLE 1 

Model properties. Tabulated from left to right are respectively stellar mass, M*, in Mq, stellar radius, R*, in Rq, 
effective temperature, t^g (k), logarithm of surface gravity, g = gmi,/r^, escape speed at the surface, 

Vescfi =Vasc(R*) = V^IGMTAR*, PHOTOSPHERIC density, Po(G CM"^), ROOT-MEAN-SQUARED AMPLITUDE OF TRANSVERSE FLUCTUATION 
AT THE PHOTOSPHERE, (Svq) (kM S""^), MINIMUM WAVE PERIOD, T^in, MAXIMUM WAVE PERIOD, Tmax, AND LOCATIONS OF FLUX TUBE 

EXPANSIONS, Ri AND i?2. 




Fig. 2. — Schematic picture of a magnetic flux tube that is 
super-radially open. The solid lines denote magnetic field lines. 
The open field region expands more rapidly than the spherical ex- 
pansion shown by the dashed lines. Our simulations treat such 
super-radially open fiux tubes. 

gation and dissipation of the waves and the consequent 
acceleration and heating of the surrounding gas, in a 
self-consistent manner. We set the outer boundaries at 
~ 25 Ri, at which the outgoing condition of both materi- 
als and MHD waves is properly prescribed (SI06). In this 
paper we neglect centrifugal force, which is a reasonable 
assumption since low mass MS stars and RGB stars gen- 
erally rotate slowly, and radiation pressure, which is also 
appropriate, where it might become important in more 
evolved asymptotic giant branch stars with dusts in the 
atmospheres. 

Radial field strength, Br, is given by the conservation 
of magnetic flux as 



BrT fir) — const. 



(1) 



where /(r) is a super-radial expansion factor. A ID flux 
tube is fixed by the field strength at the photosphere, 
and f{r). We set 5^,0 and /(r) based on our pre- 
vious solar studies (SI05; SI06). While the average field 
strength is 1 — 10 G on the solar photosphere, magnetic 
field lines are swept into inter-granulation lanes so that 



the localized strength becomes > 100 G to a few kG. 
These field lines rapidly expand and becomes < lOG in 
the coronal region. Finally in the outer region (> a few 
Ri,) the field lines distribute almost radially. 

We use the the same Brfl — 240 G and the same total 
expansion factor, /tot = 240, for all the Models. Note 
that these values give i3r,o//tot = 1 G, which is the total 
magnetic flux of the open field regions divided by the stel- 
lar surface. The average field strength becomes larger by 
adding the contributions from closed field regions ( M.ip . 
We use the same functional form of the expansion as in 
SI06 (see also Kopp & Holzcr 1976): 



/i^tot exp(l^) + h /2.tot cxp(^) + h 



cxp(r-i^i) + 1 



cxp(r-i^) + 1 



(2) 



1 



(/i 



1) exp( 



R*-Ri 



and /2 



1 



where /i 

(/2,tot — 1) exp( ^*^^^ ). A flux tube initially expands by 
a factor of /i_tot around i?i — < 7i ^ j?i + gj, c orrespond- 
ing to the 'funnel' structure (|Tu et al.l 120051 ) by closed 
loops, and followed by /2,tot times expansion around 
R2 — (J2 ^ Ri -\- 02 due to the large scale (dipole) mag- 
netic fields. The total super-radial expansion is defined 
as /tot = /i,tot/2,tot. We adopt /i,tot = 40 and /2,tot = 6 
in all the Models. R\ and i?2 of each Model are sum- 
marized in Table [TJ R\ — R^, is set to be (roughly) pro- 
portional to the scale height (cx = {GMi,/Rl)~^; 
see Equation [T51 for a more specific form) because closed 
loops are supposed to be taller in lower gravity stars if 
the surface magnetic field conditions are the same (but 
see §4.ip . Accordingly R2 is also set to be larger in more 
evolved stars. We adopt the 'widths' of the expansions, 
CTi « 0.6(i?i - Ri,) and cr2 « 0.6(i?2 - Ri,). Actually, 
the most important parameter that determines the stel- 
lar wind structure is /tot and the other parameters, 
R2, CTi, (72, /i,tot, and /2,tot, only give smaU corrections. 

We dynamically solve ideal MHD equations with ra- 
diative cooling and thermal conduction. Conservation 
of mass, momentum (both radial and tangential compo- 
nents), and energy, and evolution of magnetic field are 
respectively expressed as 



dt ^ r^fdr^ ^ 



0, 



dVr 



dp 
dr 



1 



d 



pvl 



d 



2r^f dr 



Swr^f dr 
I 2 n GM. 



(r'fBl) 



(3) 



(4) 
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dt 



e + 



47r 



87rp 



J_d_ 

j,2 J Qj, 



„2 J Qj. 



{r^fFc) + qn^O, (6) 



1 d 

ryTf dr 



(7) 



where p, v, p, B are density, velocity, pressure, and mag- 
netic field strength, respectively, and subscript r and _L 
denote radial and tangential components; ^ and ^ de- 
note Lagrangian and Eulerian derivatives, respectively; 
e — TfZi^ is specific energy and we assume the equa- 
tion of state for ideal gas with a ratio of specific heat, 
7 = 5/3; G is the gravitational constant; = kqT^^^^ 
is thermal conductive flux by Coulomb collisions, where 
Kq « 10~^ in c.g.s unit; is radiative cooling. In this 
paper we assume the solar elemental abundance for the 
radiative cooling term. In sufficiently low density re- 
gions with p < 5 X 10~^^ g cm~^ we use optically thin 
radiative loss, qn = UpUeA, where Hp and rig are pro- 
ton and electron number densities and A is the tab- 
ulated cooling function by iLandini fc Monsignori-Fossil 
(|1990[ ). In denser regions, which generally corresponds 
to low chromosphere, we adopt empirical radiative cool- 
ing, (7R = 4.5 X 10^ p (erg cm~'^s~^) derived from obser- 
vations of the solar chrom osphere (jAnderson fc Athav I 
Il989t iMorivasu et al.ll2b04l ) to take into account the op- 
tically thick effect. For simplicity's sake we switch off 
the cooling if temperature drops below 3000 K, which 
sometimes happens in the RGB stars. We might have 
to take into account effects of neutral atoms or dusts in 
such low temperature, though wc do not do explicitly in 
this paper (see ^4.3p . 

We adopt the second-order MHD-Godunov-MOCCT 
scheme (Sano & Inutsuka 2007 in preparation; see also 
SI06) to update the physical quantities. In this scheme 
each cell boundary is treated as a discontinuity, and 
for the time evolution we solve the nonlinear Riemann 
shock tube problem with magnetic pressure by using the 
Rankine-Hugoniot relations. Therefore, entropy gener- 
ation, namely heating, is automatically calculated from 
the shock jump condition. A great advantage of our code 
is that no artificial viscosity is required even for strong 
MHD shocks; numerical diffusion is suppressed to the 
minimum level for adopted numerical resolution. 

2.3. Surface Conditions 

We describe the conditions of the stellar photospheres, 
which are the inner boundaries of the simulations. 

2.3.1. Density 

The density, po, at photosphere (defined as optical 
depth « 1) is subject to stellar parameters through var- 
ious opacity effects in the atmosphere. The gas pres- 
sure, po, at photosphere has a positive dependence on 
gravity, po (x g^-^, for a given Tcs (§9 of Gray 1992). 
Pq has a negative dependence on T^s, and in the cool 
star condition (4000 K< T^s < 6000 K) the relation is 



roughly po cx T^^^^ (Figure 9.21 of Gray 1992). Then, 



Po{(x Po/Tcff) can be scaled as 



(8) 



Po (X 

10~^ g cm~^ on the sun and po's in 
the other cases are determined by this scaling relation; 
po's in Table [1] are thus given^. 



We firstly set po 



2.3.2. Wave Amphtude 

The fluctuation amplitude, Svq, at the photospheres 
can be estimated from convective flux. The generation 
of sound waves on late-type stars have been extensively 
studied (Lighthill 1952; Stein 1967; de Loore 1970; Ren- 
zini et al. 1977; Bohn 1984; Shibahashi 2005). Recently, 
numerica l experiments of s urface convection are also car- 
ried out (jStein et al.ll2004 ). These works show a similar 
qualitative trend that the excited acoustic flux, i<"a,o(erg 
cm^^s^^), has a negative dependence on surface gravity 
and a positive dependence on stellar effective temper- 
ature. In thi s pape r we adopt the relation derived by 
iRenzini et all ()1977[) . 



-Fa,o(« Po{Svl)Cs) CX g 



(9) 



where 5vq is wave amplitude, Cs is sound speed, and () 
denotes time-average. 

Strictly speaking, 5vq should be longitudinal (paral- 
lel to the field line; i.e. radial direction) component. 
In this paper we assume both transverse and longitudi- 
nal components have the same amplitude at the photo- 
sphere. This assumption is roughly satisfied at least on 
the sun because observation of solar granulations shows 
that both transverse an d longitudinal amplit udes are 
similar, w 1 - 2 km s^^ (jHolweger et al.lll978[ ). We set 
{5vq) — 1.0 km s^^ in the present sun (Model I), where 

{Svq) = \/{Svf^^). For the other stars, we can use the 



following scaling. 



{Svq[ 



fX „-0-657^7.25 



which is derived from Equations 



-,1/2. 



(10) 

and ^ (note oc 



{Svq) 's in Table[l]are determined by this relation'^ 



2.3.3. 



Wave Spectrum 

The typical period, r, of generated waves is also re- 
lated to stellar parameters, r can be scaled by eddy 
(granulation cell) turn-over timescale, ~ l/cs, where I is 
a typical granulation size, if stochastic processes of the 
convection dominate in the wave generation. Vs are dif- 
ferent in different stars, and are probably proportional 
to the pressure scale heights, Hp « c^/g. Then, we have 



T r^l/CsCC Hp/Cs 



cjg cx g-'T^-' 



off ' 



(11) 



which is consistent with recent numerical simulations 
{e.g. Stein et al. 2004). 

We use Equation pT|) to give the realistic spectra of 
fluctuations at the photosphere. In this paper, we adopt 

^ po's in Table [T] are slightly diflferent (< 10%) from the val- 
ues derived from Equation ^ This is because po's are set by the 
model atmosphere that is used to derive the approximated scaling 
of Equation JSJl. 

^ (<5i'o)'s in Table [T] are also slightly different from those based 
on Equation l|10| l. This is again because we are using the model 
atmosphere when setting Svo. 
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a broadband power spectrum of P(i^) oc ly ^ with re- 
spect to frequency, i> — l/r, where the normahzation 
is given by (Svq) ~ J^^"^^ dvP{v). For the present sun, 

we set Tmin(= = 1 Hliu and Tmax(= l/^'min) = 32 

min, giving the (logarithmic) mean value ~ 5 min, which 
dominates the solar oscillations. Twin's and Tmax's of the 
other stars can be determined by the scaling relation of 
Equation pTjl as shown in Table [TJ It should be noted 
that a shape of P{v) affects wind structure little as long 
as the same Tmin and Tmax are adopted (SI06). 

2.3.4. Wave Mode 

In this paper we only consider the open field regions 
in which the (unperturbed) magnetic fields are oriented 
to the vertical direction at the surfaces. In such a situ- 
ation, if we restrict waves traveling in the vertical direc- 
tion, the transverse components of surface fluctuations 
produce Alfven waves, while the longitudinal component 
excites compressive waves, e.g. acoustic waves, or more 
generally, magnetosonic waves. The uncompressive char- 
acter enables the Alfven waves to propagate long dis- 
tance. On the other hand, the longitudinal waves are 
quite dissipati ve because both nonlinear sh ock formation 
(ISuzuki Il2004[ ) and coUisionless dissipation (jSuzuki et al.l 
[200I " are effective. Thus, the compressive waves cannot 
travel a sufficiently long distance to contribute to the 
acceleration of the winds (Suzuki 2002). Actually, we 
tested a role of longitudinal fluctuations in several mod- 
els and found that its effect is limited in the regions near 
the surfaces (<C 2 R^,). Therefore, we only show the 
results of transverse fluctuations at the photospheres in 
this paper. 

3. RESULTS 

We initially set up the static and cool atmospheres with 
the photospheric temperatures; neither stellar winds nor 
hot coronae are set up in advance. We start the simu- 
lations by injecting the transverse fluctuations from the 
photospheres, which excite Alfven waves (see the elec- 
tronic edition for the mpeg flle of Model VI which is the 
movie version of Figure [5]) . The atmospheres are accel- 
erated as transonic winds in all the six Models by the 
dissipation of the Alfven waves f ^3.4.2p . This shows that 
the transonic outflow is a stable and natural consequence 
of the Alfven wave-driven winds of both MS and RGB 
stars. The stellar winds are settled down to the quasi- 
steady states after > twice of the sound crossing times of 
the entire simulation boxes (the sound speeds are slower 
than the Alfven speeds in our simulations) , although the 
RGB cases (Models III- VI) show large fluctuations in the 
winds by thermal instability fi i3.2p . 

First we would like to briefly mention the main results 
of our simulations. Figure [3] presents the structure of 
temperature, T (K) (top), outflow velocity, Vr (km s~^) 
(middle), and density, p (g cm^'^) (bottom) on radial 
distances from the stellar centers in unit of Rq. The 
quantities are averaged over sufficiently longer durations 
than Tmax of each Model. Figure[3]shows the typical wind 
parameters, temperature, T(K) in the atmosphere (top), 
wind speed, V2o_r» (km s~^), at 20 i?* (middle), and mass 
loss rate, M = 47r(pr^w)2o_R„ {Mq yr"-'^), evaluated at 20 
i?* (bottom), as functions of Wcsc,o; evolved stars are in 
the left side. Note that the super-radial expansions of the 
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Fig. 3. — Time-averaged stellar wind structure of the six Models. 
From the top to the bottom, temperature, T (K), radial outflow 
velocity, Vr (km s~^), and density, p(g cm"''), are plotted. The 
solid lines are the results of the 1 Mq stars; the black, blue, green, 
and red lines are the results of Models I, II, III, and IV, respectively, 
The dashed lines are the results of the 3 Mq stars; the black and 
light-blue lines are the results of Models V and VI. 

flux tubes flnish in r < 3R^, so the flux tubes are radially 
oriented ai r — 20i?*. In the top panel, we are averaging 
T{r) in 1.5 - 15 i?*. The solid error bars correspond to 
the temperature ranges of the time-averaged structure 
in this region. The dashed error bars in all the panels 
correspond to the maximum and minimum values from 
the snap-shot results f i)3.2p . 

One of the most important results is the disappearance 
of steady hot coronae in the RGB stars (Models III - VI). 
It is clearly seen in the top panel of Figure |4] that the 
average temperature drops suddenly from T ~ 7 x lO^K 
in the sub-giant star (Model II) to T < lO^K in the RGB 
stars. The large dashed error bars in the same panel 
indicate that the temperatures of the RGB atmospheres 
show large variations. One may also see in the middle 
and bottom panels that the velocities and mass loss rates 
of the RGB stars largely fluctuate. Particularly, M's of 
the RGB stars vary 2-4 orders of magnitude. Therefore, 
the red giant winds are highly structured, in contrast 
to the steady coronal winds from the MS and sub-giant 
stars. 

The middle panel of Figure [3] shows that nearly-static 
regions are formed above the photospheres in the RGB 
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Fig. 4. — Typical wind parameters of the six Models as functions 
of surface escape velocity, fesc,o (km s~^). The circles and triangles 
represent the time-averaged values of the 1 Mq (Models I - IV) 
and 3 Mq (Models V &; VI) stars, respectively. From the top 
to the bottom temperature, T (K), outflow velocity, evaluated at 
r = 20 -R* arc plotted. T (K) is averaged over the range of 1.5 -R* < 
r < 15 -R*. The solid error bars in the top panel correspond to the 
maximum and minimum values of the time-averaged structure in 
this range. The dashed error bars are the results of the maximum 
and minimum values derived from the snap-shot results ( ij3.2l l. The 
hatched region in the top panel indicates the thermally unstable 
region. 



stars, and the acceleration of the winds essentiaUy starts 
from several stellar radii. Accordingly, their wind speeds 
in the outer regions are considerably slower than the es- 
cape velocities at the surfaces (the middle panel of Fig- 
ure [31), which is consistent with the observed trend {e.g. 
Dupree). The mass loss rate rapidly increases from MS 
to RGB. The average M's of the stars with logg = 1.4 
are 10~^ — 10~ ^ M(Dyr~^, which are c omparable to the 
observed level (| Judge fc Stencellll99ll ). 

3.1. Disappearance of Steady Coronae 



We explain the mechanism of the disappearance of the 
steady corona in more detail. The main reason is that the 
sound speed (~ 150 km s~^) of « 10^ K plasma exceeds 
the escape speed, Vesdf) = yJ^GMZjr, at r > a few i?^,; 
the hot corona cannot be confined by the gravity any 
more in the atmospheres of the RGB stars. Therefore, 
the material flows out before heated up to coronal tem- 
perature. However, the temperature gap at Wosc.o — 300 
km s"^ (the top panel of Figure 2]) cannot be explained 
only by this reason. 

Thermal instability plays an important role in the 
drop of the temperature. The radiative cooling func- 
tion of the optically thin gas wit h the solar abundance 
(|Landini fc Monsignori-Fossil 1199 0) shows a decreasing 
trend on temperature in T > 10^ K, indicating that gas 
is less cooled as it is heated up. Then, thermal instability 
sets in and gas in this temperatu re range cannot stably 
exist (iKovama fc Inutsuka I [200^ . In solar and stellar 
coronal situations, however, thermal conduction plays a 
role in stabilization; drastic heating to T ~ 10^ K is 
suppressed by the downward conducti on from the up - 
per corona to the lower chromosphere (jHammer Ill982f) . 
Then, gas with T > 10^ K becomes thermally stable, 
while gas with 10^ K < T < 10^ K, which is shown by 
the hatched region in Figure [H remains unstable. 

Let us think about comparison of Cg and Vesc at a few 
i?i of the 1 Mq stars. In Models I and II, Vesc exceeds Cg 
(~ 150 km s~^) of 10® K plasma. Under such a condition 
the coronal temperatures are determined by energy bal- 
ance among heating, radiative cooling, and downward 
thermal conduction in the transition region which d i- 
vides the chromosphere and the corona (jHammer 
However, in Model III(logg = 2.4) Vesd^R*) 110 km 
s-^ and in Model IV(log5 = 1.4) Vcsc{3R*) ^ 60 km 
s 

T 

fore, the gas streams out before heated up to the stable 
temperature (> 10® K) because of the insufficient con- 
finement by the gravity. Those temperatures are ther- 
mally unstable, and then, the actual temperatures be- 
come lower than the gravity limited values. This is the 
reason why the average temperatures of Models HI — VI 
are not lO'^ — 10® K but < 10^ K. In conclusion, the ther- 
mal instability, in addition to the gravity effect, yields 
the sharp decline of the temperature from T > 10® K to 
T < 10^ K with the stellar evolution. 

The results (Figures [3] &|31) show that the temperatures 
in the winds of the 3 Mq RGB stars are systematically 
higher than those of the 1 Mq counterparts. This is 
because the energy inputs from the surfaces are larger 
in more massive stars. When fixing surface gravity, a 
star with a larger mass has higher effective temperature, 
which gives a larger photospheric amplitude (Equation 
nop . Then, more energy is injected from the surface so 
that the larger heating is achieved. 

3.2. Structured Red Giant Winds 

After the disappearance of the steady corona, the 
atmospheres and winds of the RGB stars show time- 
dependent behaviors because of the thermal instability. 
In this section, we examine the snap-shot structure of 
the red giant winds. We firstly show the results of Model 
VI (3 -M0), which is redward of the DL but more active 



^ which correspond to sound speeds of plasmas with 
6 X 10^ K and ~ 2 x 10^ K, respectively. There- 
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Fig. 5. — Snap-shot wind structure (solid) of Model VI at 
t = 6909(hr), compared with the time-averaged structure (dashed). 
From the top to the bottom, (km s~^), T (K), and p (g cm""^) 
are plotted on (r — Ri,)/ Ri,. The arrows in the middle and b ottom 
panels show the magnetized hot bubble which we inspect in i|3.2.1l 
and Figure |6] The mpeg movie which shows this simulation from 
t = to 9800 hr is available in the electronic edition. The movie 
file presents (r>x) (km s~^) in addition to (km s"'^), T (K), and 
p (g cm-3). 

than the 1 Mq counterpart (Model IV). After that, we 
discuss the other RGB cases. 

Figure O shows the snap-shot wind structure of Model 
Yl &% t — 6909 (hr) (solid) in comparison with the time- 
averaged structure (dashed) that is the same as in Figure 
[31 The mpeg movie of the simulation is also available in 
the electronic edition. Figure [5] illustrates that the sim- 
ple picture of layered atmosphere, photosphere - chromo- 
sphere - transition region - corona - wind from below, 
does not hold in the RGB star. A characteristic feature is 
that a number of hot bubbles with low densities are dis- 
tributed in cool background material. For example, in an 
inner region [r — Ri, < R^,) the two bubbles with the peak 
temperatures > 10^ K are seen; in i?^ < r—Ri, < 10 i?^, a 
couple of warm bubbles with T > 10^ K are formed. The 
hot and warm bubbles and the cool background materi- 
als are connected by the transition regions, at which the 
temperature drastically changes because of the thermal 
instability. 

The densities of the hot bubbles are lower than the 
ambient media to satisfy the pressure balance. As seen 



in the bottom panel of Figure [5] the density fluctu- 
ates typically by 2-3 orders of magnitude in the wind, 
which is related to the fact that the gas mainly con- 
sists of hot (> 10^ K) and cool (< 10"* K) compo- 
nents. Strictly speaking, however, magnetic pressure 
also needs to be taken into account for the force balance 
(t j3.2.ip . The outflow speed fluctuates as well to fulfill 
the mass conservation relation. The red giant wind is 
not a steady outward stream but an outflow consisting of 
many small-scale structures. Our simulation shows that 
both hot plasma and cool chromospheric wind coexist 
in the RGB sta r, which is consistent with observations 
of hybrid sta rs fflartm ann. Dupree. fc Raymond I Il98d : 
iHarper et an [T995: Avr es et al.lll998D . 



3.2.1. Magnettzed Hot Bubble 

How are these hot bubbles formed ? How are the high 
temperatures maintained ? We study the formation and 
evolution of bubbles by inspecting one typical hot bub- 
ble that is indicated by the arrows in Figure [5l This 
bubble is initially created at r = 1.04 Ri, dX t — 6755 
hr. It comoves outwardly with the background flow and 
eventually cooled down to « lO'* K at t ~ 7070 hr at 
r = 1.6 i?* (see the mpeg movie of Figure [5] available in 
the electronic edition). 

Figure [6] presents the zoomed-in structure just after 
the formation of the bubble (the lower panels) with 
schematic cartoon (top). The location of the shock front, 
which plays a role in the heating of the bubble, is indi- 
cated by the dotted line at r = 1.047 i?* in the panels. 
The top and middle panels showpB(= B'^/Sw; magnetic 
pressure), p, Vr^T and p. In the bottom panel, we present 
heating rate per mass. The solid line indicates the net 
heating rate that is calculated from local change of spe- 
cific heat, de/dt, and the dashed line shows the heating 
by the only MHD process which excludes radiative cool- 
ing and thermal conduction from de/dt. 

As clearly seen in the bottom panel the heating rate 
is not smoothly distributed but largely fluctuates be- 
cause the heating mainly occurs at the steepened fronts 
of waves (see i 33.4.2|) . Radiative cooling is also important; 
in some portions the net heating becomes negative (the 
solid lines disappears) in spite of the positive MHD heat- 
ing (dashed line), because the radiative cooling exceeds 
the MHD heating. This large variation of the heating 
rate by waves and cooling triggers the multi-temperature 
stellar wind. 

Pb increases at the shock front from the upstream 
(right) side to the downstream (left) side, which indi- 
cates that this is a fast MHD shock. One may find the 
density peak just behind the shock front. This is because 
an effective ratio of specific heats is small 7eff ^1-1 ow- 
ing to the thermal conduction; the density jump at the 
shock becomes large in small 7off circumstances so that 
the material should be concentrated in the narrow re- 
gion behind the sh ock to satisfy the mass conservation 
(iKim fc Rvu1[200l . As a result, the radiative cooling 
dominates the MHD heating owing to the high density, 
and this region is cooled down to form a cool dense shell 
(the bottom panel of Figure [H]) . 

Behind this cool dense shell, the hot region, which was 
heated up by the fast MHD shock, remains. The fast 
MHD shock plays an important role in maintaining this 
hot bubble. Because the magnetic field is generated by 
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Fig. 6. — Zoomed-in structure of the magnetized hot bubble 
at t = 6755(hr) (lower panels), which is the initial phase of the 
bubble indicated by the arrows in Figure |5] and schematic picture 
of a magnetized hot bubble (upper cartoon). The dotted lines in 
the lower panels show the location of the shock front. What are 
plotted in the lower three panels are as follows: top : The dashed, 
dot-dashed, and solid lines are magnetic pressure, pe = B^/8n, 
gas pressure, p, and radial velocity, Vr- The values of p and pe are 
on the right axis (dyn cm~-^) and the value of is on the left axis 
(km s~^). middle : The solid and dashed lines are temperature (K; 
left axis) and density (g cm~^; right axis), bottom : The vertical 
axis shows heating rate (erg g^^s"'^). The solid line is the net 
heating derived from the change of specific heat and the dashed 
line is the pure MHD heating which excludes the effect of radiative 
cooling and thermal conduction from the net heating. 



Fig. 7. — Mass loss rate (upper) and radial velocity (lower) of 
the stellar wind of Model VI as functions of time. Both quantities 
are evaluated at r = 20 i?* . 

the shock, the magnetic pressure in the bubble is larger 
than in the ambient gas. Therefore, the hot bubble is not 
supported by the gas pressure but by the magnetic pres- 
sure. Even though the gas pressure decreases by an oc- 
casional cooling, the bubble can avoid being compressed. 
The density is kept low to reduce the radiative loss as well 
as to increase the heating rate per mass. The hot bubble 
survives much longer time than ambient (average) cool- 
ing time until the magnetic topology changes; the high 
temperature is maintained around 315hr, while a typical 
cooling time, TcooI, of gas with T — lO^K (A w 10~^^erg 
gj^3g-i. q2.2p is 



7*0001 
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1 min 



10 
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(12) 



To summarize, hot bubbles are initially created by the 
variations of heating rate, and their high temperature is 
maintained thanks to the support by magnetic pressure. 

An intriguing concept of magnetized hot bubbles is 
that hot plasma can be heated up even in open field re- 
gions. It has been generally considered that closed loops 
are necessary to confine hot plasma in the atmosphere of 
a RGB star since the gravity is not sufficiently strong. 
However, our simulations show that the hot bubble can 
exist rather long time (^^ 10 days) by magnetic topol- 
ogy, i.e. a fast MHD shock, in open field regions. It is 
supposed that such hot bubbles are ubiquitously created 
and they are possible sources of continuous soft X-ray 
emissions in hybrid stars as discussed later ( §3.2.3|) . 

3.2.2. Time-Variation of Wind Parameters 

Because the red giant wind consists of many small- 
scale structures, the mass loss rate and the outflow speed 
in the outer region also show time-dependent behaviors. 
Figure [7] presents M {Mq yr~^) (upper) and V20 (km 
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s""'^) (lower). The upper panel shows M varies more than 
three orders of magnitude from 1 x 10^^*^ to 5 x 10"^ 
Mq yr~^. The lower panel also shows that V20 largely 
varies from 10 km s^^ to 60 km s~^. 

The jumps are seen in M and W20fl* a.t t = 7200, 8750, 
and 9020 hr. These jumps indicate that dense blobs pass 
at r = 20 i?* at those times. We examine one of the 
blobs which passes at t — 7200 hr (This blob is seen in 
Figure [5] as well). This is a shock wave originally formed 
in an inner region. The shock amplitude has been ini- 
tially small, but it grows as the shock moves outward 
by sweeping-up a number of smaller dense blobs as well 
as by amplification in the stratified atmosphere with the 
decreasing density. Finally the compression ratio of the 
shock increases to more than 10. (This value exceeds the 
adiabatic limit = 4 since our simulations, which take 
into account cooling and thermal conduction, are not 
adiabatic.) Before this shock arrives, M is almost con- 
stant and V20Rj, shows smoother feature. After the shock 
reaches 20i?^, both M and W20fl, show fluctuating be- 
haviors. This is because the downstream region contains 
larger perturbation than the upstream region. 

3.2.3. Radiative Flux 

Hot bubbles in RGB atmospheres become sources of 
soft X-ray and extreme UV (EUV) radiation that is ob- 
served in hybrid stars. We present the radiative fluxes, 
normalized by the bolometric luminosity, from the entire 
star of Models VI as functions of time in Figure [S] The 
blue, green and red lines are the radiative fluxes from 
the gas with T < 2 x lO'' K (optical; chromosphere), 
2 X 10'* K < T < 5 X 10^ K (UV; transition region) and 
5 X 10^ K < T (EUV/soft X-ray; corona). We caU them 
cool, warm, and hot components, respectively. We also 
plot the radiative flux of the gas with T > 10® K, which 
is a part of the hot component and corresponds to soft 
X-ray emission. 

When deriving the radiative fluxes from the entire star, 
we assume spherical symmetry: 

i(Ti <T<T2)=A^ j drr^qnin <T< T^), (13) 

where Ti and T2 are the minimum and maximum tem- 
peratures of each component. While our simulations only 
treat the open magnetic field regions, in reality a certain 
fraction of the surface is covered by closed structures. 
Equation (jl3p assumes that the radiative fluxes from 
both open and closed regions have the same r depen- 
dence. Although this treatment is not strictly correct, 
we can give rough estimates. 

According to the upper panel of Figure [3 the hot 
component shows strong intermittent activities, while 
the cool component, which dominates the total radiative 
flux, stays roughly constant within 10~^ < Lcooi/^* < 
10^*. The behavior of the warm component is be- 
tween the cool and hot ones; when the hot compo- 
nent is switched-off or stays in low activity phases {e.g. 
7500-8000 hr) Lwarm/i* ~ 10"'', but once the hot com- 
ponent becomes active, iwarm/i* increases to > 10~^ by 
the contributions from the transition regions of hot bub- 
bles. 

It is also shown that the hot bubbles with T > 10® 
K exist during the only half of the presented duration. 




5500 6000 6500 7000 7500 8000 8500 
t(hr) 

Fig. 8. — Upper. Normalized radiative fluxes, L/L-t, of cool 
(T < 2 X 10* K; blue), warm (2 X 10* < T < 5 X 10^ K; green), and 
hot (T > 5 X 10^ K; red) components of Model VI as functions of 
time elapsed from the beginning of the simulation. The radiative 
flux of T > lO'' K gas, which is a part of the hot component, is 
also plotted (dashed line). Lower: Radiative flux of the soft X-ray 
(T > 10® K) that takes into account the absorption by the cool 
wind materials (black solid), in comparison with the original flux 
without the absorption (pink dashed). 

However, we should note that this is based on the spheri- 
cal symmetry approximation. In reality, hot bubbles are 
created in different open flux tubes at different times, 
and they are ubiquitously distributed in the atmosphere 
to be continuous soft X-ray sources. 

When discussing observed X-ray flux, we need to take 
into account the absorpt ion by Hydro gen in the outer 
cool chromospheric win d (|Avresll200l . Then, the ob- 
served flux of the hot component would be smaller than 
that shown in the upper panel of Figure [S] In particular, 
soft X-ray radiation is supposed to severely suffer the ab- 
sorption since most of the hot bubbles with T > 10®K 
are distributed in the inner region (r < 2 i?^). 

In the lower panel of Figure[8]we compare the observed 
soft X-ray flux (T > 10® K) that takes into account the 
absorption (solid) with the original flux (dashed). Ten- 
fold obscuration of the soft X-ray emissions sets in at co l- 
umn density, TVh = 5 x 102^-5 x lO^* cm'^ (jAvres I2005D . 
We adopt the most optimistic case of TVh = 5 x 10^^ 
cm""^: 

/drr^^RlO-^H/Sxio-cm-^ (^4) 
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Fig. 9. — Temperature structure (upper panels) and time- 
variations of radiative fluxes (lower panels) of Models III (left pan- 
els) and IV (right panels). The snap-shot results of temperature 
(solid lines) are at t = 1295 and 6050hr in Models III and IV, re- 
spectively. The time-averaged results that are the same as in Figure 
[3] are also shown by the dashed lines. The three components of the 
radiation is the same as in Figure |8] 

When deriving A^h from the simulation result, we only 
consider the gas with T < 8000 K to pick up neutral 
Hydrogen and assume the solar abundance. The lower 
panel of Figure [8] illustrates that most of the soft X- 
ray emissions are absorbed. Around t ~ 6500 and 8000 
hr, the soft X-ray could be observed, and the Lx/Li, is 
reduced to < 10^^ from 10"^ — 10~*. Interestingly, this 
value, Lx/Li, < 10~^, is typical for hybrid stars (e.g. 
Ayres 2005), although our calculation, which is based on 
ID and do not consider detailed radiative processes, is 
too simplified for further quantitative arguments. 

3.2.4. Other RGB Cases 

We have examined the time-dependent behaviors of the 
3 Mq red giant wind (Model VI). In this subsection we 
study the 1 Mq stars on both blueward (Model III) and 
redward (Model IV) sides of the DL. Figure [5] shows the 
snap-shot temperature structure and the time-variation 
of the radiative fluxes of Models III and IV. The up- 
per panels (temperature) show that the atmospheres and 
winds of the IMq RGB stars also consist of many hot 
and warm bubbles embedded in the cool chromospheric 
materials. It is confirmed that RGB winds are generally 
structured with many blobs and bubbles. 

The temperature of bubbles are lower in the more 
evolved star (Model IV) . This trend is also clearly seen in 
the plots of the radiative fluxes (the lower panels). The 
emission from the hot gas (T > 5 x 10^ K) exists only 
during the short term (t = 5800 hr) in Model VI, while 
the hot bubbles with T > 5 x 10^ K almost always ex- 
ist in Model III. Namely, once a star with 1 Mq crosses 
the DL, EUV/soft X-ray emissions from hot plasma with 



r > 5 X 10^ K happen to disappear. However, the quali- 
tative wind features, hot/warm bubbles in cool materials, 
are similar in Models III and IV. Both temperature and 
the number of hot bubbles decrease with the evolution 
from Models HI to IV, because the heating rate per mass 
decreases on account of the increase of the density in the 
atmosphere (Figure 

The emission levels of the hot components of the 1 Mq 
stars are smaller than that of the 3 Mq star. This is be- 
cause Svq (Table [1]) of the surface convection is smaller 
in a less massive star owing to the lower Tcff (Equation 
[Tn|) . Our result is consistent with the observed trend 
that UV/X-ray emission o f RGB stars are positively cor- 
related with stellar mass (iHiinsch fc Schrdder |[l996l) . 

The average emission level of the hot component of 
Model III is Lhot/-^* ~ 10~^, which is marginally smaller 
than the typical observed X-ray flux {Lx/L-i, — 10^^ — 
10~^; e.g. Ayres 2005) from RGB stars blueward of the 
DL. Moreover, the emissions are mostly from the bub- 
bles with T < 10®(K) in our simulation, corresponding 
to EUV rather than soft-X ray emissions, and the ef- 
fect of the absorption will further reduce the observed 
flux. Then, it would be difficult to explain the observed 
X-ray from RGB stars before crossing the DL by only 
the hot bubbles in open field regions. Probably, hot 
plasma, which is efhciently confined in closed magnetic 
loops, also contribu te to the X-ray emission. Incidently, 
iRosner et al.l (|1995f ) speculated that the DL was reflec- 
tion of the change of magnetic topology; as a star crosses 
the DL, the magnetic structure becomes open and the 
hot plasma is liberated to be X-ray deficient (see H4.ip . 

3.3. Slow Red Giant Winds 

The wind velocity also drastically decreases along with 
the stellar evolution as shown in Figure [H As briefly 
stated previously, this is mainly because the nearly static 
regions are formed above the photospheres in the RGB 
stars and the winds are essentially accelerated from the 
effective surfaces located several We see this in more 
detail here. In order to investigate the drastic change of 
the wind speed, we compare Models near the transition 
at fosco = 200 — 300 km s^^ (the top panel of Figure^]). 
Figure [10] shows the simulated density structure (solid) 
of Models II (left), V (middle), and III (right), in com- 
parison with the density structure based on hydrostatic 
equilibrium (dotted) and that derived from assumed con- 
stant velocities (dashed). 

Hydrostatic equilibrium gives 



p = pi exp 



= Pi exp 



V "cff 



1 1 

r Ti 



91 n 



-*off 



(15) 



where ri is a certain reference position, z = r — ri, 
gi = GMi,/rf, and pi is the density at r = ri. We 
here assume the "isothermal" atmosphere with an effec- 
tive sound speed, UeS = + B'j^/STrp, which takes into 
account magnetic pressure associated with Alfven waves 
in addition to gas pressure, where a is isothermal sound 
speed. For Models II, V, and HI, we adopt — 74, 55, 
and 43 km s""'^, respectively, which are estimated from 
the simulation results. We set ri at the coronal base 
(Model II) or the locations where the average temper- 
atures start to increase (Models HI & V). Note that 
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Fig. 10. — Density structure of Model II (left), V (middle), and 
III (right). The results are the same as in Figure |3] but compared 
with the densities derived from constant wind speeds (dashed lines) 
and from the hydrostatic equilibrium (dotted lines) (see text). 

ri 2:; ( 1.01 - 1.02)i?^ in these Models (see Figure [TT] in 
t j3.4.ip . Figure [TU] shows that the densities in the sub- 
sonic (inner) regions can be approximated by the hydro- 
static structure, which is consistent with the result of 
general stellar wind theory {e.g. section 3 of Lamers & 
Cassinelli 1999). This indicates that the atmospheres are 
supported by both magnetic and gas pressure. 
The dashed lines plot 

P^{r'f)'\ (16) 

which is derived from substituting constant Vr into the 
steady-state continuity equation of mass, pVrT^ f =const. 
The starting point is set in the low chromosphere in each 
Model. Because in Model II (the left panel) the density 
decreases faster due to the stronger gravity, the dashed 
line is far above the solid as well as dotted lines. This 
indicates the acceleration of the wind should start near 
the surface. 

On the other hand, in Models III (the right panel) 
the decrease of the density is slower since the gravity is 
weaker. As a result, the solid line follows the dashed 
line in the inner region, which indicates that the wind 
is not accelerated, but stays almost a constant speed. 
The solid line is gradually diverging from the dashed line 
around i?eff ~ several i?*, which we call an "effective" 
surface, and at which the wind starts to be practically 
accelerated. In r < i?efi, the average wind speed is kept 
much slower than 10 km s~^(^ Cg), and the atmosphere 
is almost static, as seen in Figure [31 the acceleration is 
not effective there because of too much material. The 
longitudinal motions seen in this region in the top panel 
of Figure [3] are not by outflow but due to compressive 
waves (! j3.4.2p . The wind structure of the other RGB 
stars (Models IV and VI) is essentially similar to that of 
Model III. 

The gravity is smaller at r = i?off than at the surface, 
hence the slow wind can escape outwardly. Therefore, 
we can expect that the terminal velocity is not an order 
of the surface escape speed but an order of the escape 
speed, Wosc(^cfi), at r = RcS- This can explain why 
observed wind speeds are much smaller than the surface 
escape velocities in RGB stars {e.g. Dupree 1986). In 
Model III we can set i?cff — 8i?^, hence, Wosc(^Gff) — 70 
km s~^. This value is an order of the terminal velocity 
which is supposed to be slightly faster than the obtained 
= 40—50 km s^^ when the effect of the acceleration 
in r > 2QRi, is taken into account. 



The situation of Model V (the middle panel) is be- 
tween Models II and III. The simulated density follows 
the trend of the non-acceleration at the beginning. How- 
ever, the solid line departs from the dashed line earlier 
than in Model III. In this case the wind becomes grad- 
ually accelerated from a few R-^. Setting i?off ~ 3 i?*, 
we have Wesc ~ 150 km s~^ which is an order of the ob- 
tained V2{) Kt. ■ The difference between Models III and V 
arises from the difference of the radii; although these two 
Models have the same g, the decrease of the density is 
faster in Model V when measured in r/i?^,, {z = r — Ri, 
is larger in Equation [TS]) . 

In concluding, we have the natural conclusion that the 
wind speed is determined by the escape speed at r = i?efi 
where the wind practically starts to flow out. R^g does 
not coincide with the stellar surface but moves outward 
as a star expands. Therefore, the wind speeds of the 
RGB stars are slower than the surface escape speeds. 

3.4. Energetics 

The main energy that drives the stellar winds is the 
turbulent energy of the surface convection in our simula- 
tions. The Alfven waves that are excited by the footpoint 
motions of open flux tubes carry the energy outwardly 
and the dissipation of the waves leads to the heating 
and acceleration of the stellar winds. In this subsection, 
we quantitatively investigate the energetics of the stellar 
winds driven by Alfven waves, and finally determine the 
mass loss rate in a forward manner. 

3.4.1. Heating 

We study the structure of lower atmosphere because it 
is important not only in determining the basal properties 
of the stellar winds but in controling the reflection of the 
Alfven waves. In Figure [11] we show the time-averaged 
density (top) and temperature (middle) which are the 
same as in Figure[3|but plotted on (r — i?*)/i?i to enlarge 
the inner regions. Note that the physical distance {e.g. 
km) of the horizontal axis becomes larger for a star with 
larger i?^. In the bottom panel plasma (3 values defined 
as 



are plotted. 

In the MS (Model I) and sub-giant (Model II) cases, the 
steady corona lies above the (geometrically) thin chro- 
mosphere and these two regions are divided by the sharp 
transition region. In the RGB cases (Models III- VI), 
the temperature gradually increases outwardly from the 
chromosphere. This increase of the temperature is differ- 
ent from the formation of "warm corona" but the result 
of the time-average of intermittent creation of hot bub- 
bles; the location at which the temperature increases is 
the starting point of the formation of hot bubbles. How- 
ever, it is still worth studying the average trend from an 
energetics point of view. 

The lower panel illustrates that the variations of /3's are 
qualitatively similar in the stars with logg = 2.4 — 4.4; 
the P values decrease outwardly to ^ 10~^ but turn to 
increase from certain locations. This is because when the 
magnetic pressure dominates the gas pressure by ^1000 
times, the dissipation of a small fraction of magnetic en- 
ergy gives an enormous effect on the gas. Indeed, the 
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Fig. 11. — Structure of the inner regions of the six Models. From 
the top to the bottom, T (K), p (g cm~'^), and plasma /3 values are 
shown on (r — iJt)/i?*. The results in the top and middle panels 
are the same as in Figure [3] The solid lines are the results of the 
1 Mq stars; the black, blue, green, and red lines are the results 
of Models I, II, III, and IV, respectively. The dashed lines are the 
results of the 3 Mq stars; the black and light-blue lines are the 
results of Models V and VI. 



energy transfer from the magnetic fields to the plasmas 
is carried by the dissipation of the Alfven waves in our 
simulations and the gas is heated up before /3 becomes 
further smaller. As the temperatures rise, the decreases 
of the densities become slower (Equation [15]) and /3's in- 
crease due to the decreases of i?'s. The densities at which 
the average temperatures start to increase are controlled 
in order that the plasma /3 values do not become too 
small. The geometrical depth of the inner chromosphere 
(the region between the photosphere and the location of 
temperature rise) is larger in a lower log 5 star since the 
decrease of density is slower. 

In the stars with \ogg = 1.4, the /3 values do not be- 
come as small as the other cases since the decreases of 
the densities in the low chromospheres are too slow to 
achieve such small /3's. Therefore, strong heating do not 
occur and the average temperatures are lower than the 
higher gravity stars. 

In the 3Mq stars the temperatures start to increas 
from more inner positions than in the IMq counterparts. 
This is because a more massive star gives larger pho- 
tospheric amplitude and input energy of Alfven waves 
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Fig. 12. — Adiabatic constants, S^'s (erg cm~^s~^), of outgo- 
ing Alfven (solid), incoming Alfven (dotted), and outgoing slow 
(dashed) MHD waves in Model III. 

(Equation [TU]) . Then, hot bubbles are heated up from a 
deeper [i.e. denser) region. 

3.4.2. Dissipation and Reflection of Waves 

In SI05 and SI06, we claimed that the nonlinear gener- 
ation of MHD slow ( « sound) waves is the main dissipa- 
tion process in our ID simulation of the solar wind; the 
energy of Alfven waves is transferred to the slow waves, 
which dissipate by shocks after steepening of the wave 
fronts. We show that the same mechanism operates in 
RGB stars at least under the ID approximation. We 
firstly inspect the dissipation processes by using the re- 
sults of Model III, and after that, we compare the four 
cases of the 1 Mq stars (Models I- IV). 

We are injecting the broadband spectra of the surface 
fluctuations and in the upper regions not only outgoing 
but incoming waves of different modes travel simultane- 
ously. Hence, the actual simulations show turbulence- like 
feature {e.g. see v± in the movie of Figure [5] available 
in elec tric edition). A techri ic of the mode decompo- 
sition (iCho fc Lazarianl[2003l) is a useful tool to study 
such turbulent phenomena. Figure [T^] plots the following 
adiab atic constants derived from wave action (jJacaues I 

[ilz3), 

fph rif{rc) 

for outgoing Alfven, incoming Alfven, and outgoing slow 
MHD (sound) waves of Model HI as functions of dis- 
tance, where Sv^^ and Vph are amplitude and phase 
speed of each wave mode. In our simple ID geometry, 
Alfven waves are purely transverse and slow MHD waves 
are purely longitudinal; the energy flux of Alfven waves 
is calculated from v±^ and i3_L, and the energy flux of 
slow MHD waves is from fluctuating components of Vj. 
and p. A characteristic feature of incoming (outgoing) 
Alfven waves is that v±_ and B±^ are in (in opposite) 
phase. The decomposition of each wave mode can be 
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done by using these facts^. The subscript, c, denotes 
that S is normalized at rc where the flux tubes expand 
by fc = 40, to exclude the effect of the adiabatic loss. 

Sc can be used as a measure of the dissipation. Al- 
though Sc and energy flux are identical in static media, 
they are different in moving media. Sc is conserved in 
expanding atmosphere if a wave does not dissipate, while 
it is not the case for the energy flux. For the incoming 
Alfven wave, we plot the opposite sign of Sc so that it 
becomes positive in the sub-Alfvenic region. 

Figure [12] shows that the outgoing Alfven waves dissi- 
pate quite effectively; Sc{= \Sc\) becomes only ~ 10^^ 
of the initial value at the outer boundary. First, a sizable 
amount is reflected back quite near the surface, which is 
clearly illustrated as the incoming Alfven wave that is 
following the outgoing component with a slightly smaller 
level. This is because the wave shape is considerably 
deformed owing to the steep density gradient. In this re- 
gion, the temperature is small, T ~ 3000 K, which gives 
a short pressure scale height, Hp (Figure [TT|) . Then, the 
wavelength of the Alfven waves becomes comparable or 
longer than Hp, which results in the effective reflection 
of the waves oore "eTanflOQl . In this Model 111, more 
than 60% of the initial energy flux of the Alfven waves 
are reflected back before reaching the upper warm region 
with the average temperature, ~ 10^ K (^ ^3.4.3^ . 

Second, slow MHD waves are generated in the corona 
(see Sakurai et al.2002 for solar observation) as shown in 
Figure [121 The generation of the slow waves is triggered 
by the variation of the magnetic pressure^ Bi/87r, with 
the Alfven waves (Kudoh & Shibata 1991). This ex- 
cites density perturbations, which are slow MHD waves 
in magnetically dominated plasma. Slow waves are gen- 
erally longitudinal, and then, they eventually su ffer non- 
linear steepening to lead to shock dissipation (I Suzuki I 

[2002) . The density fluctuations of slow waves also work 
as mirrors to Alfven waves. The variation of the density 
directly leads to the variation of the Alfven speeds (ex 
B/^/p), which enhances the reflection of the Alfven waves 
by the deformation of the shape. These reflected waves 
further nonlinearly interact with the pre-existing outgo- 
ing waves, which plays a role in the dissipation. 

The process discussed here is the parametric decay in- 
stability of Alfven waves due to three-wave (outgoing 
Alfven, incoming Alfven, and outgoing slow waves) in- 
teractions (Goldstein 1978; Terasawa et al. 1986). The 
decay instability is not generally efficient in the homo- 
geneous background since it is a nonlinear mechanism. 
Therefore, it has not been considered as a major pro- 
cess of the acceleration of the solar and stellar winds. 
However, the density gradient of the background plasma 
totally changes the situation. The amplitude, v±, of 
Alfven waves is inevitably amplified with the decreas- 
ing density to conserve S. On the other hand, the 
Alfven speed, 'yA(c>c B/^/p), does not increase so much 
because the decrease of Br (flux tube expansion) com- 
pensates the decrease of p. As a result, the nonlinearity 

This method is strictly valid only for linear and short wave- 
length waves of which amplitudes are sufficiently smaller than 
phase speed and wavelengths are sufficiently shorter than a vari- 
ation scale of background phase speed (i.e. WKB approximation 
holds). For nonlinear and long wavelength waves, which we want to 
deal with, this treatment is only statistically true llCho fc LazariaiTl 

[2003 ) . 
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Fig. f 3. — Properties of Alfven waves in the 1 Mq stars (Models 
I - fV). The top panel exhibits the adiabatic constant, Sc, and the 
bottom panel shows the nonlinearity, (fi)/fA- The black, blue, 
green, and red lines are the results of Models I, If, fit, and fV, 
respectively, 

((u^) /wa) increases as the Alfven waves travel upwardly 
(see Figure [T5)) . which leads to efficient dissipation. 

Figure [T2I compares the dissipation of the Alfven waves 
in the four Models of the 1 Mq stars. In the upper panel 
S'c's are compared, and in the lower panel nonlinearities, 
{v±}/va, are shown, where {v±) = \J (u^). The figure 
shows that the wave energies drastically attenuate just 
above the surfaces in higher gravity stars because of the 
reflection in the chromospheres. In these stars the den- 
sities decrease more rapidly in the inner chromospheres, 
leading to the increases of wa's. Then, the wavelength of 
the Alfven wave with the typical frequency (e.g. several 
min. in Model 1) becomes longer than Hp 1/g which 
itself is substantially smaller for larger g. This enhances 
the reflection as explained above. 

Furthermore, in stronger gravity stars, the chromo- 
spheres are geometrically thin (r — i?^, <C R*)- Other- 
wise if the chromosphere extends to r ~ R^, > R^, as in 
the RGB stars, a curvature effect increases Hp by r/R^, 
(Equation [T5|l. which reduces the wave reflection. This 
is also a reason why the reflection becomes relatively im- 
portant in the higher gravity stars. Finally, only less 
than half of the input wave energies can go beyond the 
inner chromospheres in Models 1 - 111 ( ^3.4.3|) . 

Let us turn to the waves in the outer regions. Wave 
dissipation is more rapid in lower gravity stars. This 
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Fig. 14. — The fractions of wind kinetic energy, radiative loss, 
and wave energy reflected back to the surfaces of the 1 Mq stars 
converted from the input energy of Alfven wave as functions of 
Vesc,o (km s-i). 

is because the densities in the atmopsheres are larger 
in those stars due to the longer Hp and accordingly 
the Alfven speeds, cx 1/y^, becomes slower. In such 
a condition Alfven waves more easily become nonliear 
{{v±)/vA ^ 1 in the bottom panel of Figure [T3|) . which 
enhances the wave dissipation. Then, one may expect 
that the wave energies are more effectively transferred to 
the winds in more evolved stars. However, the situation 
is not so simple because in these stars radiative losses 
become efficient due to the higher densities. We study 
the wind energetics to discuss these issues in the next 
subsection. 

Before closing this subsetion we would like to remark 
on dissipation lenghts of the Alfven waves. Since the dis- 
sipation processes of Alfven waves are quite complicated, 
a simpliphied phenomcnological approach has been often 
adopted to treat Alfven waves in solar and stellar wind 
models {e.g. Hartmann & MacGregor 1980); an action 
of Alfven waves is assumed to follow 

5cxexp(^-^-y^^ , (19) 

where I is a (constant) dissipation length {— 0.1 — 10 i?*). 
We can directly estimate / from our simulations. The 
upper panel of Figure [13] illustrates that the decreases 
of ScS are roughly on straight lines in the log(r/i?*) — 
log Sc diagram. This means that the dissipation lengths 
increase on r. For example, in Model IV Z < 1 i?* in 
r < 3 i?*, but / ~ 10 Rn in r > 10 the assumption of 
a constant dissipation length is not appropriate. 

3.4.3. Energy Transfer from Alfven Waves 

A fraction of the initial wave energy is used for the ac- 
celeration of the stellar winds, while the rest of them is 
lost by radiation and returned back to the stellar surface 
by counter propagating waves. We estimate kinetic en- 
ergy fluxes of the stellar winds at the outer boundaries, 
integrated radiative losses, and wave energy fluxes leak- 
ing inward from the stellar surfaces, respectively, in unit 



of erg/s : 

{pVr^fr^Ut, (20) 

(•'•out 

/ drr^fqR, (21) 

{p{Svl_)vA).Rl (22) 

where Sv±_- denotes amplitude of incoming 
Alfven waves, and subscripts and 'out' indicate 
that the values are evaluated at the stellar surfaces and 
the outer boundaries, respectively. Here, we do not 
discuss energy carried by waves, enthalpy, and thermal 
conduction away from the outer boundaries because 
their contributions are smaller. Figure fT4l presents these 
energies normalized by the input energy fluxes from the 
surfaces, 

{p{Svl)vA).Rl (23) 

as functions of Uesc,o- 

As one can see, the fractions of the wave energies that 
are converted to the wind kinetic energies are small. Typ- 
ically, only ~ 1 % of the surface turbulent energies are 
transferred to the stellar winds. However, these tiny frac- 
tions are enough to keep the observed M and Voo ~ W20fl* 
of the winds from MS and RGB stars. 

Most of the energy is lost by the radiation and the 
wave reflection. In the less evolved stars (Models Till), 
the energy losses are dominated by the wave reflection as 
discussed in ii3.4.21 the densities decrease more rapidly 
because of the higher gravity, and then, the wavelengths 
become longer than the scale heights, which results in 
the effective reflection via the deformation of the wave 
shape. The rapid decreases of the densities also make 
the radiative cooling less effective; it is important only in 
the chromospheres (T < lO'^ K) and low coronae where 
the densities are sufficiently high (volumetric radiative 
cooling is cx in optically thin and cx p in optically 
thick gas). This also leads to the relative dominance of 
the wave reflection in the total energy loss. On the other 
hand, the trend of the most evolved star (Model IV) is 
opposite since the decrease of the density is more grad- 
ual; radiative cooling becomes relatively important first 
because it is efficient even in the outer region and second 
because the Alfven waves do not suffer the refiection so 
much. 

Let us turn back to the kinetic energy part. The ten- 
dency on Wosc,o is not monochromatic. From Model I 
(MS; Wcsc,o = 617 km s^^) to II (sub-giant; Wcsc.o — 347 
km s~^), the fraction of the kinetic energy increases 
mainly because the mass flux, pvr^ increases. However, 
the fraction drops almost by one order of magnitude from 
Models II to HI (wcsc.o — 195 km s^^). This is because 
of the drastic decrease of the wind speed (Figure |4]). 
Although the density increases from Models II to HI, 
the effect of the velocity dominates since the kinetic en- 
ergy flux is proportional to pv^. From Models HI to IV 
(^csc.o = 110 kni s^^), the decrease of is small, and 
the fraction of the kinetic energy again increases simply 
by the increase of the density. 

3.4.4. Mass Loss Rate 
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Fig. 15. — Comparison of the mass loss rates of the simulations 
(vertical axis) with the Reimers' relation (filled symbols; horizontal 
axis) and Equation II31I I derived from the scaling relations of wind 
energetics (open symbols; horizontal axis). The circles and the tri- 
angles indicate the results of the IA/q and SAf© stars, respectively. 

The bottom panel of Figure [4] shows the mass loss rate 
of the IMq star rapidly increases by ~ 5 orders of mag- 
nitude from MS to RGB phases. However, the increase 
of M is continuous, which is in contrast to the velocity 
and temperature showing the gaps between the sub-giant 
and RGB stars. This is because the density in the wind 
drastically rises during this phase, which compensate the 
rapid decrease of the wind speed. 

Since our simulations adopt a forward approach, M 
is not an input parameter but an output directly deter- 
mined from the simulations. Therefore, it is straightfor- 
ward to study how M evolves with the stellar evolution 
from the wind energetics. We begin by comp a ring o ur 
results with the empirical relation bv iReimers I (|1975( ): 



J W Reimers — 4 X lU T^R , , , , , , . Mq 



[AU/Mq) 



/yr, 



(24) 

where 77^ is a proportional coefficient of an order of unity. 
This is based on a dimensional analysis of stellar wind 
energetics. The filled symbols in Figure [15] compares 
the obtained M with the MRcimers for ?7r = 1. The 
comparison indicates that the Reimers' relation roughly 
gives reasonable estimates. However, smaller 77r(< 0.1) 
seems to be systematically favored in the higher gravity 
stars. To search a better relation, we investigate the 
mechanism that controls M from the energetics of our 
Alfven-wave driven wind simulations. 

In the wave-driven winds, the radiation part {L-^) 
should be replaced by wave energy in Equation ([M]) . In 
an open flux tube, a fraction, e, of the initial energy flux 
of input Alfven wave is transferred to the stellar wind 
where e is exactly what has been derived in iJ3.4.3l (Fig- 
ure [H]). This can be written as 

,,2 



'fpvr- 



(25) 



where i^w,star = Pq{Svq)va is wave energy flux. The left- 
hand side is evaluated at the surface and the right-hand 



side is at a sufficiently distant location at which / = 
/tot and V w Voc- Multiplied by An, a transformation of 
Equation (|25p gives 

„2 



eLw = M- 



(26) 



where — 47ri?jFw//tot is wave luminosity and M — 
Anpvoor'^- Here Fw//tot is used instead of because 
the fraction of the area occupied by open field regions at 
the surface is l//tot- Then, an equation that determines 
M becomes 



2e- 



2e- 



;,0 



:,0 -^w 

(27) 



where we add subscript 'w' to distinguish from AfRcimers 
and the simulated M. Although e varies on Vcsc.o, we 
have almost constant 



^^IcJ^lo ~ 0.02, 



(28) 



because f^/^esc positively correlated with e as shown 
in Figure [HI 

As for the wave luminosity, we can use the scaling re- 
lations of Equations ^ and (fTU)) . and then. 



i^w//tot (Brfi/ftot) VPoi^'^o) 



1.4 X 10 



CIT 

(ff/ffo) 



(29) 



where a is the Stephan Boltzmann constant, subscript, 
0, denotes the solar values, and we are using the same 
Br.o/ftoti= 1 G) in aU the Models. 

Substituting Equations (pS)) and into Equation 
([?f|) . we finally have 



1.7 X 10-^% 



(L./L0)(i?./i?0) (Teff/Teff,0)S 



(30) 



1.7 X lO-'V 



{LjLQ){RjRQf{T,s/T,s,er 



(31) 



(M,/M0)2 

in M0yr~^, where we leave a degree of freedom in r]^{^ 
1), similarly to the Reimers' relation (Equation [M)) since 
e is not a constant. Note that the dependences of g and 
Toff in Equation ([50]) are from the depende nces of the en- 
ergy fl ux of the generated sound waves by iRenzini et al.l 
(Il977l) (Equ at ion [9]); if we used other results by e.g. 
IStein I (l967D . iBohnI (|1984D . which give slightly differ- 
ent dependences, the scaling relations of Af^ would be 
also slightly different. Afw does not have dependence 
on the properties of magnetic flelds since we assume a 
constant -Br,o//tot(= 1 G) in Equation (p9|) . In more 
realistic situations, however, i?r,o//tot would show star- 
to-star variation, and Af^ would have dependence on the 
magnetic fields. 

The open symbols in Figure [15] compares the simulated 
M and M^. One can see that most of the data points 
are distributed near rj^ = 1. While the Reimers' formula, 
mainly targeting more luminous red giants, does not give 
good Af estimates for the less evolved stars (filled sym- 
bol; Equation IM)) . our scaling gives a better explanation 
to these stars. Then, the scaling of Equation ((3T|l can 
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connect the mass loss rate of less evolved stars to that of 
more evolved ones that can be handled by the Reimers' 

f ormula. 

ISchroder fc: Cuntzl ()2005D also proposed a formula of 
mass loss rates of wave-drive winds, focusing on more 
evolved stars. In their formula, additional dependences 
on g and T^s also appear, while they are different from 
ours. 

4. DISCUSSION 

4.1. Magnetic Fields 

When performing the simulations, we have to set up in 
advance the surface fluctuations and the magnetic flux 
tubes. While we can reasonably estimate the ampli- 
tude and period (spectrum) of the surface fluctuations 
from the conditions of the convection zones ( !j2.3[) . there 
are large uncertainties in the properties of the magnetic 
fields. 

We only consider the open magnetic field regions 
because our aim is to study the stellar winds. We 
assume the same magnetic field strength (-Br,o = 
240 G) at the surface and expansion factor (/tot = 
240) in the different stars. Although it is unclear 
how in reality the open magnetic flux tubes are, we 
had better consider effects of different flux tu be s (see 
iFalceta-Conyalves, Vidotto. fc Janenco-Pereira I ()2006D 
for configurations of flux tubes) . 

We have actually studied the effects of the strength and 
configuration of magnetic flux tubes on the solar wind by 
both numerica l simulations (SI06) and analytical method 
(|Suzukill200"6[ ). We showed that Br,o//tot IS a more im- 
portant parameter in controling solar wind properties 
than individual i?r,o and /tot, if nonlinear Alfven waves 
dominantly work in the accelerati on. This is also con - 
sistent with recent observation by iKojima et al.l (|2005( ) . 
For example, flux tubes with different B^.o and /tot but 
the same ratio -Br,o//tot give similar wind properties if 
the photospheric amplitudes are the same. An input en- 
ergy flux is ~ p{5vq)'^va,o o: Brfi at the surface, and it 
is diluted according to l//r^ in an expanding fiux tube. 
The acceleration of the stellar wind mainly occurs in the 
outer region where the super-radial expansion of the fiux 
tube already finishes. Therefore, i?r,o//tot, which is pro- 
portional to Sc («energy fiux) of the input Alfven wave 
measured at the outer region, becomes an important con- 
trol parameter in the Alfven wave-driven wind. More- 
over, a mass flux, i.e. M, is determined almost solely by 
{Svo)'^, and very weakly dependent on i?r,o//tot (SI06). 
Then, i?r,o//tot plays a role only in tuning the terminal 
velocity of the solar wind. 

This argument is directly applicable to red giant winds. 
M's of the stellar winds (the top panel of Figure 2]) are 
determined with little uncertainty since they are mainly 
determined by {Svq)'^ of the surface convections. How- 
ever, when we discuss the wind terminal speeds, i?r,o//tot 
needs to be taken carefully; wind speed becomes faster 
in a fiux tube with larger i?r,o//tot- 

Sr,o//tot(= 240/240 = 1 G) corresponds to the total 
magnetic fiux over the whole open field regions divided 
by the entire stellar surface. To discuss the average field 
strength we have to take into account closed regions, 
which occupy a sizable fraction of the surface. If 10% 
of the surface is occupied by open field regions, the aver- 



age strength becomes 10 times, namely lOG. This value 
is quite reasonable from an energetics point of view. The 
magnetic energy density at the surface is B^q/Stt ~ 4 
erg cm~'^, which is considerably smaller than the kinetic 
energy of the surface turbulence, p{Svq) ~ 10'^ — 10"* erg 
cm~^, of the simulated stars. It is sufficient only if ~ 0.1 
% of the turbulent energy is converted to the magnetic 
energy by motions of charged particles {i.e. current) in 
the surface convec tion. 

Related to this, iRosner et all ()199lL Il995f ) tried to in- 
terpret the DL by the change of magnetic topologies. 
The stellar rotation becomes slow with the evolution so 
that the generation of magnetic fields by large-scale dy- 
namo is switched-off. Magnetic fields are generated only 
by (small-scale) turbulent convection in the surface lay- 
ers in evolved stars. As a result, closed fields, which 
have dominantly covered the surface in unevolved stars, 
become open in evolved stars redward of the DL. There- 
fore, hot plasma, which has been well confined by closed 
loops, would fiow out. 

If this scenario is correct*^, surface field strength, i?r,o, 
becomes smaller in evolved stars. However, we expect 
that the field strength in the outer regions (oc -Br.o//tot) 
is less affected, because an decrease of Br a firstly reduces 
a fraction of closed structures rather than strength of 
open fields in the outer region; even if B^fl decreases. 
Bo I f tot stays more or less constant before all the field 
lines become open. Because i?r,o//tot is more important 
than Brfl as stated above, we anticipate that our results 
of the stellar winds are affected little even if the magnetic 
field topologies change with the stellar evolution. 

4.2. Previous Wind Calculations 

iHartmann fc MacGregor I (|1980l ) investigated stellar 
winds from late-type stars with explicitly taking into ac- 
count the effects of Alfven waves. They calculated several 
stars with different wave energies and dissipation param- 
eters under the steady-state condition from assumed in- 
ner boundaries. They concluded that the dissipation of 
Alfven waves resulted in coronal heating in higher grav- 
ity stars (MS and sub-giant) and mass loss in lower grav- 
ity stars (RGB). Our analysis of the time-average wind 
structures basically confirm this tendency, whereas our 
simulations have revealed the various important time- 
dependent aspects ( g3.2p . 

Most of previous works claim that it is difficult to ac- 
celerate slow and dense red giant winds by Alfven wave s 
(jHolzer et al. II1983I : ICharbonneau fc MacGre"gMl ll995)^ 
The main problem is that a sizable amount of the in- 
put wave energy remains in the supersonic region, which 
leads to faster wind. There is also a similar problem 
in the concep t of a supersonic transition locus, which 
iMuUan I (|1978f) introduced to explain the onset of massive 
wind. The main difference of our results from the previ- 
ous works is that the large subsonic regions are formed 
in our simulations due to the nearly static regions above 

^ It should be noted that lAvres, Brown, fc Harperl I I2003I) re- 
cently proposed a different scenario for the DL: even though closed 
loops, which are sources of X-rays, are distributed in stars redward 
of the DL, these X-rays are obscured by larger column density of 
above chromospheric material. This density effect is a reasonable 
explanation of the DL. 

^ Although these works focus on more massive RGB stars than 
ours, the result will not change if less massive stars are considered. 
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the surfaces. Therefore, most of the input Alfven wave 
energies are deposited in the subsonic regions, which re- 
sults in the slow winds, although the dissipation is not 
much faster than in the previous model calculations. 

We speculate that our dynamical wind s olutions corre- 
spond to "inner solutions" introduced by iHolzer et al. I 
(|1983f ). They pointed out that there are two transonic 
solutions in some cases; one is the inner solution in which 
the wind is accelerated very gradually and the wind speed 
is slower, although they mainly studied the other one, 
the outer solution, in which the acceleration starts just 
above the surface. Our dynamical simulations automat- 
ically select stable transonic wind solutions. Then, our 
results might show that the inner solution is the stable 
branch of Alfven wave-driven winds. 

We would like to remark on an issue of time-dependent 
winds. Our simulations show that some portions of the 
RGB atmospheres transiently accrete {vr < 0), which 
lasts longer than the assumed longest wave periods {e.g. 
see the dashed line in the top panel of Figure [5|) . This 
phenomenon may be an example of the hysterisis transi- 
tion betwee n subs o nic ac c retion and supersonic wind dis- 
covered by IVelli I ()1994[ ): iDel Zanna. Velh. fc Londrillo I 
(11998) • 

4.3. Limitations 

Here we discuss processes and effects which need to be 
studied in more detailed calculations. 

4.3.1. Heating of Low Chromosphere: Sound Waves? 

We do not input longitudinal fluctuations at the photo- 
spheres, which would directly drive sound waves, because 
they do not contribute to the heating and acceleration of 
the stellar winds {e.g. Judge & Carpenter 1998) (In our 
simulations sound waves are generated from Alfven waves 
in the upper regions ; ^3.4.2p . However, the sound waves 
directly from the surface might be important in the heat- 
ing of low chromosphere (e.g. Carlsson & Stein 1997). 
We switch off the radiative cooling in the inner regions if 
T < 3000 K and p > 5 X lO^i^ g cm~^ by the technical 
reason fi i2.2p . By taking into account the longitudinal 
perturbations at the photospheres, we may carry out a 
more self-consistent treatment of the low chromospheres 
without switching off the radiative cooling artificially. 

4.3.2. Multt- dimensional Effects - Wave Generation- 

In our simulations we only consider the wave gener- 
ation at the photospheres in the ID flux tubes. How- 
ever, interactions between flux tubes may be important 
in generation of waves, though 2D or 3D modeling is 
required to study such processes. For example, recon- 
nection events between closed loops could be sources of 
small flare-like events, and these events probably excite 
compressive wa ves at locations above the photosphere 
(lSturrocklll999l) . 

The magnetic reconnection between a closed loop and 
an open field line also possibly excites A lfven waves trav- 
eling upwardly (|Mivagoshi et al.l l2004f ). These waves 
would also play a role in the acceleration of red giant 
winds. Further detailed analysis is necessary to estimate 
how such waves produced by activities involving closed 
loops are important with respect to the energetics. 



4.3.3. Multi- dimensional Effects - Wave Propagation- 

We treat the propagation of Alfven waves in the ID 
flux tubes. In other words, Alfven waves are assumed to 
travel in the same way in neighboring flux tubes. How- 
ever, if Alfven speeds are slightly different in neighboring 
field lines, the waves become out of phase even if they 
are generated in the same way at the surface. Then, 
these Alfven waves diss ipate by phase mixing th rough 
resistivity and viscosity (jHevvaerts fc Prie"stlll983f l. 

Alfven waves in stellar winds become more or less 
turbulent-like. Then, Alfvenic turbulence cascades to 
higher frequency in the transverse direction with respect 
to u nderlying magnetic fields (iGoldr eich fc Sridhar I 
I1995D . Therefore, Alfven waves might dissipate by dif- 
ferent channels in more realistic 3D calculations, while 
in our ID simulations the Alfven waves dissipate mainly 
through the decay instability (i j3.4.2p . 

4.3.4. Collisionless Effects 

We assume the MHD approximation, namely a mean 
free path with respect to Coulomb collisions, Zmfp = 

9.38 X 10' cm ^ylgg cm-^ ' smaller than a typical scale. 
Here, we can take the minimum wavelength, Amin — 
UATmin, of the Alfven waves as a typical scale, hence, 

is the requirement for the MHD condition. Otherwise 
we have to consider collisionless effects of plasma. This 
condition breaks in r > 2 i?^ of the MS star (Model I) 
and in > 5 i?^ of the sub-giant star (Model II), while 
in the RGB stars (Models III - VI) it holds in the entire 
simulation regions because the densities are high and the 
temperatures are low. 

In the collisionless regime, compressive (slow/fast 
MHD) wave s suffer collisionless damping (Barnes 196^ 
ISuzuki et al.l 120061 ). Then, the dissipation rate of the 
slow MHD waves, which are nonlinearly generated from 
the Alfven waves, may be changed. This might further 
modify the mode conversion rate from the Alfven waves. 
Then, the dissipation of the Alfven waves might be indi- 
rectly affected in Models I and II. 

4.3.5. Atoms, Molecules, and Dusts 

In the RGB atmospheres, the temperatures sometimes 
become 10'^ — 10^ K as shown in the snap-shot results 
(Figures[5]&[9]). Although even in such low temperatures 
a fraction of the gas is expected to be ionized by the 
UV/X-ray emissions from hot bubbles, a sizable amount 
of neutral atoms remain, or even molecules can also exist 
if T < 3000 K (Tsuji 1973). The snap-shot resuhs also 
show that the temperatures partly decreases further to ~ 
1000 K, and in such circumstances dusts are also formed. 

The existence of atoms and molecules indicate that the 
gas is not fully ionized but partially ionized. The MHD 
approximation implicitly assumes that the coupling be- 
tween ions and neutrals is sufficient so that the gas can 
be treated as one fluid. Otherwise, the friction between 
these two components plays a role in the dissipation of 
Alfven waves. We can check the condition of ion-neutral 
coupling by c omparing wave frequency and a ion-n eutral 
collision rate ([Draine. Roberge. fc Dalga rno |[T983[ ). We 
have found that the coupling is sufficient with respect 
to the low-frequency waves we are dealing with; we do 
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not have to worry about the effect of weakly ionized gas. 
However, if high frequency Alfven waves are generated 
by turbulent cascade, these waves dissipate by the ion- 
neut ral damping more effectively than in fully-ionized 
gas (Holz er et aLl ll983l. 

Molecules can be used as diagnostics to understand 
physical conditions of RGB star atmospheres. Observa- 
tions of CO show that quasi-static cool (T ~ 1000 — 2000 
K) regions containing molecules are in the atmospheres of 
RGB stars, typically at ~ 2i?* (jTsuii II1988I: iTsuii et all 
|1997() . Our simulations of the RGB stars also show the 
quasi-static regions in the atmospheres and the temper- 
atures partly becomes low enough to form molecules. 
The cooling by mo lecules drive thermal instabihty a t 
T = 2000 - 3500 K (|Muchmore. Nuth. fc SteiiceTlfigsl . 
similarly to the ions in T > 10^ K. This will be impor- 
tant in controlling thermal properties of the atmospheres 
and winds of RGB stars. 

Radiation pressure on dus ts is not efficie nt in 

driving winds of RGB stars (lJudge fc Stencel I 119911 : 
IVidotto fc Janenco-Pereiral I2006D . although it could be 
important in asymptotic giant branch stars. However, 
Alfven waves may d issipate by the cyclotron resonanc e 
with charged dusts (jVidotto fc Janenco-Pereiral I2006D . 
which might also give a correction to the dissipation of 
the Alfven waves. To study this effect, we have to care- 
fully consider the formation processes of dusts in inter- 
mittently existing cool (T ~ 1000 K) gas. 

5. SUMMARY 

We have studied the winds of the intermediate and low 
mass RGB stars near the DL, by comparing with the 
coronal winds from the MS and SB stars. These stars all 
have the surface convective layers, which play a major 
role in driving the stellar winds as well as the UV/X-ray 
activities because other mechanisms such as centrifugal 
force and radiation pressure are insufficient. The surface 
turbulences excite various types of waves, among which 
the Alfven wave travels a longer distance and contribute 
to not only the atmospheric heating near the surface but 
the wind acceleration in > several Ri,. 

We have, for the first time, successfully carried out the 
time-dependent MHD simulations of the red giant winds 
in the open magnetic field regions from the photospheres 
to the sufficiently distant locations 25 i?*). We inject 
the fluctuations from the photospheres, of which the am- 
plitudes and spectra are estimated from the properties 
of the surface convection. We have determined the ve- 
locities, densities, and temperatures of the stellar winds 
in a self-consistent manner by taking into account the 
radiative cooling and thermal conduction; the propaga- 
tion and dissipation of the MHD waves are automatically 
treated, and the subsequent heating and acceleration of 
the gas are also self-consistently handled by the jump 
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conditions of the MHD shocks, using the nonlinear MHD 
Godunov-MOCCT algorithm. 

The stellar winds are driven by the nonlinear dissi- 
pation of the Alfven waves through the conversion to 
compressive waves, whereas other dissipation processes 
might become important when performing 3D model- 
ing/simulations. The energy fraction that is transferred 
to the stellar wind against the initial wave energy is quite 
small, ~ 1%, and the rest of the energy is lost by the wave 
reflection (less evolved stars) and radiation (more evolved 
stars). However, this tiny fraction is enough to explain 
the observed stellar winds from MS and RGB stars. 

When the stars evolve to slightly blueward positions 
of the DL, the steady hot coronae with temperature, 
T > 10^ K, disappear. In the atmospheres and winds 
of the RGB stars magnetized hot bubbles and cool chro- 
mospheric gas co-exist. The hot bubbles are created in 
fast MHD shocks and supported by the magnetic pres- 
sure. Therefore, the densities of the bubbles can be kept 
low. The radiative cooling is reduced to balance with the 
heating so that the bubbles survive much longer than the 
cooling timescale. 

The red giant winds are structured with many bubbles 
and blobs, and the mass loss rate and wind speed show 
large time- variations. The hot bubbles possibly become 
sources of EUV and soft X-ray emissions of hybrid stars, 
whereas the cool chromospheric wind might absorb the 
emissions. The temperatures of the hot bubbles become 
lower in more evolved and less massive stars, which is 
consistent with the trends inferred from the observations 
of hybrid stars. 

The wind velocities also drop to 10-100 km s~^ in the 
RGB stars, which are much slower than the escape speeds 
at the stellar surfaces. This is mainly because the static 
atmospheres are formed above the photospheres and the 
stellar winds are practically accelerated from several stel- 
lar radii. The terminal velocities are regulated by the 
escape speeds there, and the slower winds can escape. 

We have finally derived a relation (Equation [3T|) that 
determines the mass loss rates of MS to RGB stars from 
the energetics of the simulated wave-drive winds. This 
has a different dependence on the stellar parameters from 
the well known Reimers' formula which mainly focuses on 
more evolved giants. The derived relation gives a better 
explanation of M of less luminous stars and can be used 
complimentarily to the Reimers' formula. 
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